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The spectral functions of the pseudospin correlation functions in the bilayer quantum Hall system 
at V = 1 are investigated numerically, where the pseudospin describes the layer degrees of freedom. 
In the pseudospin-ferromagnetic phase, the lowest-energy excitation branch is closely connected with 
the ground state through the fluctuations of pseudospin Sy and Sz, and it plays a significant role 
on the tunneling properties in this system. For the system with very small tunneling amplitude 
and layer separation smaller than the critical one, the system-size dependence of calculated spectral 
function Ay^ suggests the superfluidity on the tunneling current in the absence of impurities. 



Recently the bilayer quantum Hall (QH) system |Q] has 
often been investigated from the viewpoint of quantum 
ferromagnetism. In fact, the pseudospin can be intro- 
duced to this system by a simple definition of assigning 
the upper/lower layers to the pseudospin t / i and 
the pseudospin ferromagnetism can be realized for the 
total Landau level filling v = \ jm (m: odd integers) . In 
particular, the v = \ bilayer QH system has extensively 
been studied both experimentally and theoretically ||l|. 

The pseudospin rotational symmetry in this system is 
degraded because of the difference between the intralayer 
and interlayer Coulomb interaction for a finite layer sep- 
aration d (and of the interlayer tunneling). Since the 
two layers tend to have equal numbers of electrons to re- 
duce the interaction energy, the z component of the total 
pseudospin (half the difference in the numbers) tends to 
vanish, (^z) = 0, in the ground state. Thus the pseu- 
dospin will lie in the xy plane in the pseudospin space, 
where the system maintains an invariance under rotations 
about the z axis, i.e., the SU(2) symmetry is reduced to 

u(i) i- 

The electron correlation alone pushes the ground state 
for V = 1/to towards a ferromagnetic one so the bi- 
layer V = \/m QH system behaves like an easy-plane 
itinerant-electron ferromagnet. The U(l) symmetry is 
further degraded in the presence of the interlayer tunnel- 
ing, since the tunneling amplitude behaves like a mag- 
netic field acting on the pseudospin, as seen in the tun- 
neling Hamiltonian, = — AgAs Sx- Here AgAs is the 
energy difference between the symmetric and antisym- 
metric combinations of isolated layer states and is pro- 
portional to the interlayer-tunneling amplitude 

So a bilayer QH system is characterized by two di- 
mensionless parameters, d/l and AsAs/(e^/eZ) [I: the 
magnetic length, e: dielectric constant of the host mate- 
rial). For V = 1, the phase diagram against these two pa- 
rameters is obtained experimentally Q , and the diagram 
shows that the QH state disappears for d > dc and that 
the critical separation dc increases as AgAS increases. 

Within this QH effect region, the ground state for i/ = 
1 evolves continuously from tunncling-dominated (single- 
particle like) to correlation-dominated (many-body like) 
as AgAg is decreased, and there is no intervening non- 



ferromagnetic region between the two regimes 0. More- 
over, within this pseudospin-ferromagnetic region, the 
ground-state wavefunction is analytically given as a spin- 
squeezed state [H , where the fluctuation ASz is supressed 
by the electron correlation. 

By using a sample with very weak interlayer tunnel- 
ing (AgAs/(e^/eO ^ 10~^), Spielman et al have recently 
showed the huge resonant enhancement of the zero bias 
tunneling conductance The possibility of the exci- 
tonic superfluidity was also discussed because of the tun- 
neling spectroscopy reminiscent of the Josephson effect, 
of which controvertial predictions have been done [Q. 
Very recently, further tunneling-spectroscopy measure- 
ments were performed with the in-plane magnetic field 
applied to the system, and the observation of a linearly- 
dispersing Goldstone mode in a bilayer QH ferromagnet 
was claimed [|o[. This in-plane field effect has been dis- 
cussed in some theoretical papers JTT|-p^ . 

In case of very weak interlayer tunneling as realized 
in these experiments the U(l) symmetry of the 

system is expected to recover almost completely. Some 
simple model systems with this symmetry have been in- 
vestigated and the tunneling properties have been dis- 
cussed The Chern-S imons Landau-Ginzburg ap- 
proaches for bilayer QH systems p^-p^ were also per- 
formed for the case of vanishing interlayer tunneling. 

On the other hand, some discussions based a BCS-like 
wavefunction with a phase variable (j), which determines 
the orientation of the pseudospin magnetic moment, have 
also been done Q. This wavefunction, however, is com- 
patible with the SU(2) symmetry, and does not reflect the 
U(l) symmetry of the system as well as the Hartree-Fock 
approximation does not. Thus an approach that fully re- 
flects the system's symmetry and takes small (but non- 
zero) tunneling amplitudes into consideration has yet to 
come for the study of the tunneling-current properties in 
this system. That is exactly the purpose of the present 
paper. 

By the exact diagonalization method for small-size sys- 
tems, we obtain the spectral functions of pseudospin cor- 
relation functions numerically. We show that in the 
pseudospin-ferromagnetic phase the lowest-energy exci- 
tation branch is closely connected with the ground state 
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through the fluctuations of pseudospin Sy and Sz, and 
that it plays a significant role for the tunneling properties 
in this system. For the system with very small tunneling 
amplitude and layer separation smaller than the critical 
one, the system-size dependence of calculated spectral 
function Ayz suggests the superfluidity on the tunneling 
current in the absence of impurities. 

Our numerical calculations were done for bilayer spher- 
ical systems Q, where each surface of the spheres is 
passed through by 25" flux quanta. For the case of 
ly — 1 we consider in this paper, 2S ~ N — 1 for iV- 
electron systems. In terms of annihilation and creation 
operators for the m-th spatial orbit in the lowest Lan- 
dau level with pseudospin a, the total pseudospin op- 
erators are given as Sz = I]m("LT"™T ^ 

Sx = Em(aLT«™i + arnia™T)/2, Sy = Em(aLT«'"i - 
<^mi'^rn^)/2i. Pscudospiu Sx and Sz are related with 
the interlayer-tunneling and interlayer-voltage part of 
the Hamiltonian, respectively, as Ht = — Asas Sx and 
Hy = eVSz- Pseudospin Sy is related with the inter- 
layer current / as / oc as Sy. Our calculations were 
performed for the cases of zero interlayer voltage {V — 0). 
The spin density operators are given as 

Sxin) = + *|(0)*t(0)]/2, = 

- ^'|(r2)«'|(r2)]/2i, where a unit vector 
O denotes a position on each sphere. Their 'Fourier 
transforms' are given by = / dftYLAii^) Sa{^), 

where a = x, y, z, and Yi^f,j{^l) is a spherical harmonics 
for angular momentum L and its z-component M. In 
particular, the components for L = M = are related 
with the pseudospin operators as Sqq = Sa/V^- 

The spectral functions [ p7|jT^ ] of pseudospin correlation 
functions at zero temperature are given as 

i 

x6{uj-E, + Eo), (1) 

where |$o) and are the ground state and the i-th 
eigenstate of the Hamiltonian, respectively. The diagonal 
element AaaiL,Lu) among these spectral functions is real 
and non- negative because of their definitions in Eq. (|^). 
Among the off-diagonal elements, Ayz{L,uj) is closely re- 
lated with the finite-wavevector conductance g on the 
interlayer current because g(L,uj) oc i^^ Asas ^yziL,uj) 
from the Kubo formula. We note that Ayz{L, uj) is purely 
imaginary and that the spectral functions in Eq. are 
independent of M because of the rotational symmetry of 
the system. 

Table || shows the fluctuations of pseudospin Sz and 
expectation values of pseudospin Sx in the 1^ = 1 ground 
state for N = 10, d/l = 1.0, and t = AsAs/(eVeO = 
10"'', lO"'"^, 10~^. For such small Asas, the values are 
proportional to AsaSj because AgAS is proportional to 
the interlayer tunneling amplitude and behaves like a 
magnetic field along the x axis in the pseudospin space. 



Table | indicates that the U(l) symmetry in the easy- 
plane ferromagnet is recovered almost completely in case 
of very small Asas- In fact, for t < 10~'^, intra- and 
inter-layer correlations between electrons in the ground 
state are found to be little dependent on the value of t. 
On the other hand, the interlayer-tunneling features are 
strongly dependent on it. 

In Fig. 0, the excitation spectra and spectral functions 
in the i/ — 1 bilayer QH system are shown for N — 10, 
d/l = 1.0, and t = 10~^. The excitation spectra are 
shown by open circles in Fig. 0(a). In Fig. |l|(b) and (c), 
the strength for each excitation in the spectral functions 
Ayy {L,uj) and Azz {L,U!) is shown against the angular 
momentum L and excitation energy w by the area of the 
corresponding circle, respectively. We note that the en- 
ergy is measured from the ground-state energy and that 
the unit is /el. 

The lowest-energy excitation branch can be seen for 
< L < 3 in Fig. |l|(a). We made sure that the en- 
ergy values in this branch (in particular, the one for 
L = 0) are quite dependent on the system size N and 
that these ones decrease as N increases. This branch 
is expected to be seen as a linearly-dispersing excitation 
mode ||l^ for sufficiently large N and very small Asas- 
This pseudospin-wave excitation branch is closely con- 
nected with the ground state through the fluctuations of 
pseudospin Sy and Sz- 

As seen in Table |, the U(l) symmetry is recovered 
almost completely for t < 10^"*. In fact, we made sure 
that the result for the spectral function Axx is almost the 
same as that for Ayy shown in Fig. [^(b), and that each 
pair of excited states corresponding to Axx and Ayy, re- 
spectively, is nearly degenerate. On the other hand, as 
seen in Fig. |l|(c), most of the spectral weights for Azz are 
occupied by the lowest-energy excitation branch. This is 
also the case for Ayz, because Ayz and Azz include the 
same factor for each term in Eq. (Q) . 

In Table ||, the strength of i~^Ayz (L, u) in the lowest- 
energy excitation branch in the v — I pseudospin ferro- 
magnet is shown for = 10, d/l — 1.0, and t — 10^*, 
10^^, 10^^. It indicates that these values of strength are 
proportional to t, i.e., to Asas- Thus the tunneling con- 
ductance is considered to be proportional to the square 
of Asas, because g{L,uj) cx z~-'AsAS^iyz {L,u}). It is also 
found that the strength increases as angular momentum 
L decreases. That is, the long-wavelength part of the 
pseudospin-wave excitation mode plays a significant role 
on the interlayer-tunneling properties in this system. 

^^"^ Ayz {L,u}) /t in the 



In Table III 



the strength of i 
lowest-energy excitation branch in the v = 1 pseudospin 
ferromagnet is shown for d/l = 1.0, L = 0, 1, 2, and 
A^ = 6, 8, 10. It has a system-size dependence varying 
approximately as A^'^ (A ~ 3). This is also the case for 
d/l = 0.5. For larger layer separations such as d/l = 1.5 
and 2.0, however, such a system-size dependence is not 
obtained. In fact, for d/l — 1.5 and 2.0, the strength in- 
creases or decreases as A increases, and the system-size 
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dependence is much weaker than that for d/l = 0.5 and 
1.0. 

For sufhciently small Asas ( ^ O.Ole^/e^ ), it is known 
that the pseudospin ferromagnetism is destroyed for layer 
separation d larger than the critical one dc — 1.3Z |p|,po|. 
Among the values of layer separation adopted by us, 
d/l = 0.5 and 1.0 belong to the case of d < dc, and 
d/l — 1.5 and 2.0 satisfy d > dc- As mentioned above, 
our results for d < dc show a system-size dependence of 
the spectral function Ayz varying approximately as 
(A ~ 3), which is much stronger than usual iV- linear de- 
pendence. On the other hand, the possible superfluidity 
on the tunneling current in this system is considered to 
show an iV-square dependence. It is not clear whether 
this discrepancy between our result and predicted one re- 
sults from small system sizes in our calculations or not. 
At least, however, the system-size dependence of Ay^ in 
our results is much stronger than A'^-linear one and so 
suggests the superfluidity on the tunneling current in the 
absence of impurities. For clarifying interlayer tunneling 
properties in this system, further studies of the impurity 
effects and numerical calculations for larger systems 
are desired. 

For d > dc, such a strong system-size dependence does 
not appear in our numerical results. Significant differ- 
ences between the two cases of d > dc and d < dc are ex- 
pected to become clearer in larger systems because of the 
differences in the system-size dependence. These quali- 
tative differences between the two cases are consistent 
with the results in a recent experiment Lastly we 
note that the spectral weights of i~'^Ayz occupied by the 
lowest-energy excitation branch decrease rapidly as the 
layer separation increases. 
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TABLE I. The fluctuations of pseudospin Sz and expec- 
tation values of pseudospin Sx in the v = 1 ground state are 
shown for N = 10, d/l = 1.0, and t = AsAs/{e^/d) = lO"*, 
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For such small Asas, the values are proportional 
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16.02 X 10"* 


16.02 X lO-"* 


16.02 X lO-*" 
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10.56 X 10-* 
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8.328 X 10-" 
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4.806 X 10-* 


4.806 X lO-"* 


4.806 X 10-*= 



TABLE II. The strength of the spectral function 
Ayz {L,uj) in the lowest-energy excitation branch in the 
1/ = 1 pseudospin ferromagnet is shown for A'' = 10, d/l = 1.0, 
and t = AsAs/(eVeO = W'^, 10~^ , W'^ . The strength in- 
creases as angular momentum L decreases, and it is propor- 
tional to t. Most of the spectral weights are occupied by this 
pseudospin-wave excitation branch. 
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1.617 
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8.329 



TABLE III. The strength of r'^Ay^, {L,uj)/t in the low- 
est-energy excitation branch in the u — 1 pseudospin ferro- 
magnet is shown for d/l = 1.0, L = 0, 1, 2, and A = 6, 8, 
10, where t = Asas/ [e^ / d) . It has a system-size dependence 
varying approximately as A"^ (A ~ 3). 
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FIG. 1. The excitation spectra and spectral functions in the v = 1 bilayer QH system are shown for TV = 10, d/l = 1.0, and 
t = AsAs/(e^/£0 = W~^- (a) The energy levels are shown against the total angular momentum by open circles. The strength 
for each excitation in the spectral functions (b) Ayy (L, ui) and (c) Azz {L, uj) is shown against the angular momentum L and 
excitation energy ui by the area of the corresponding circle. The energy is measured from the ground-state energy and the unit 
is e^/el. 
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